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We consider the scenario of large scale structure formation with tau neutrino with mass in the 
MeV range and lifetime of order of years which decays into a massive majoron with KeV mass. 
The latter are the present-day cold dark matter particles. In contrast to the usual collisionless dark 
matter, the majorons have a relatively strong self-interaction and the picture of structure formation 
is rather different from the standard. The decay J — * 77 leads to the existence of an X-ray line 
at E-y = mj/2 which could be detectable. Electron and muon neutrinos are expected to be very 
light, as required in order to account for the solar neutrino deficit through v e to oscillations. 
Supersymmetry with spontaneously broken R parity provides a natural particle physics model for 
our scenario. 
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The measurement by COBE of the large angle anisotropy of the cosmic microwave background ra- 
diation has determined the magnitude of primordial density fluctuations at very large scales. This has ruled 
out the simple cold dark matter (CDM) scenario of structure formation 0, since with the power on large 
scales fixed as measured by COBE and with the assumption of a scale free primordial fluctuations spectrum 
the model leads to approximately twice more power at galactic and cluster scales than observed. If one wants 
to describe structure formation in a model with a flat primordial fluctuation spectrum one has to assume that 
the onset of matter dominance (MD) took place later than in a simple CDM model. This can be achieved e.g. 
in an open universe with h 2 Q ~ 0.2, where h = H/ (lOOkm/ sec/Mpc) and fl — pj p c . This model is disfavored 
by inflationary scenarios which (at least in simple versions) predict Q = 1. Another possibility is the mixed 
(hot+cold) dark matter scenario [|j with CIcdm — 0.7 and Qhdm — 0.3. One can also get a flat universe with 
Qtot = 1 and low matter density if the cosmological constant A is nonzero. Such models also give a satisfactory 
description of the observed structure || . 

The common shortcoming of these models is that they all demand a certain amount of fine tuning. For 
example, one would generally expect that the contribution from hot and cold dark matter to differ by many 
orders of magnitude. 

Recently there appeared a renewed interest to the idea of structure formation with unstable particles 
The recent models assume that there exists a massive long-lived particle, usually the tau neutrino with mass 
in the MeV range, decaying into massless species at the epoch when the mass density of the parent particles 
dominates the energy density of the universe. Correspondingly, the present-day energy density of relativistic 
particles is bigger than in the standard scenario and the onset of the MD stage takes place later. All these 
models have so far been based on the assumption of two unrelated components. A possible exception might be 
the models with non-zero A if a consistent adjustment mechanism is found which, as argued in ref. leads 
generically to an noncompensated amount of vacuum energy of order m 2 Pl /t 2 . 

In this letter we will consider a model for the unstable particle scenario in which both the unstable 
particles and the particles of cold (or possibly warm) dark matter are closely inter-connected. In fact, the decay 
of the former produces particles of the present-day dark matter. 

A necessary background model of this kind in particle physics was proposed some time ago in an attempt 
to find a phenomenologically acceptable way to spontaneously break R parity in supersymmetric extensions of 
the standard SU{2) ® U(l) model 0. The breaking of R parity occurs around the electroweak scale and implies 
the violation of lepton number Moreover, in the simplest realization adopted in ref. |7) lepton number 
is violated only through R parity violating interactions. As a result, the magnitude of the R parity violating 
effects is related to the large Majorana type mass of the tau neutrino, which could easily reach the MeV range. 
As shown in ref. || the v T is unstable and decays into the light plus a majoron J JToj ] . A possible breaking 
of global symmetries by gravity would give a nonzero mass to the corresponding (pseudo) goldstone boson Jll| . 
In this context it has been suggested that the majoron may pick up a mass of the order KeV, which is in the 
right range for it to play the role of cold dark matter particle [ fl2"| . 

We assume that the values of mass and lifetime of v T are such that there existed a period when it 
dominated the cosmological energy density. We denote the corresponding red-shift when it started as Z\. 
Another interesting moment was when the tau neutrino decayed and the universe returned to an RD stage. 
The corresponding red-shift is Zd- And at last there is the moment, which is especially important for structure 
formation, when the contemporary MD stage started. The onset of this stage was at red-shift Z2- The red-shift 



z\ is determined by the equality of energy densities of relativistic, p re i, and nonrelativistic, p m , matter: 

zi + 1 ~ 2.9 x 10 8 r v m M (l + O^ 4 )" 1 (1) 

Here tum — m„ T /MeV, r u is the ratio of v T number density to that of normal massless neutrinos, and £ = Xj/T 7 . 
For light v T (m < 3 MeV) i « 1 and for heavy v T (m > 10 MeV) £ ~ (4/ll) 1/3 . The value of the photon 
temperature at this moment is T\ ~ 69KeVmMr I/ (l + 0.34£ ) . This result is valid if T\ is larger than the 
majoron mass mj. We assume here that this is true. The other two interesting red-shifts are given by 

z d + 1 ~ 1.3 x 10 5 (m M r,)- 1/3 (r y )- 2 / 3 (2) 

z 2 + 1 ~ 2.8 x 10W(1 + 10 3 (m M ^) 4/3 T y 2 / 3 )- 1 (3) 

The upper bound on the universe age does not permit neutrinos to have mass above a few tens eV 
or, in the case of nonzero cosmological constant, a few hundreds eV Jl3[ | . This bound can be avoided only if 
the v T is unstable or has annihilation cross section considerably larger than in the standard model. In fact 
in the scenario we consider here both ingredients are present, so that neutrino masses in the MeV region are 
cosmologically allowed. The successful prediction for the nucleosynthesis of light elements forbids Majorana tau 
neutrino masses in the range 0.5 < hum < 35 if the latter has the standard weak annihilation cross section and 
is stable on the nucleosynthesis time scale i.e r > 100 sec |l4|, In our case, however, this bound can be 
avoided because v T has a rather strong diagonal Yukawa coupling to majorons, 

£^ } = }Jv?a e v T + {.\. (4) 

In many majoron models g — m u /V where V is the scale at which spontaneous lepton number violation occurs 
which is, by assumption, the weak scale, so that g — 10~ 5 itim / (V/100 GeV). In order to suppress r u sufficiently 
so that the massive v T would not distort the nucleosynthesis predictions on the one hand, and sufficiently large 
for a noticeable delay of the onset of the last MD stage on the other, we need g slightly above 10~ 4 . This value 
fits nicely with the assumption that the v r mass lies in the MeV range. It is also surprisingly close to the values 
assumed in the spontaneously broken R parity model Jj^j. For such g values v T v T annihilations into majorons 
substantially reduces the v T number density during nucleosynthesis and after [*]. 

Another useful restriction follows from the bound on the total energy density in the universe. For any 
form of matter we have an upper bound: p x < ptot = 10.5r2/i 2 KeV/cm 3 . Applying this bound to the massive 
majorons gives rj < O.lflh 2 /rriK where mx is the majoron mass in KeV and rj = nj/ng, where nj is the 
majoron number density and no is the number density of normal massless neutrinos, uq = 0.18T 3 . The majorons 
in our model are produced in the early universe either thermally or possibly by the lepton number violation 
phase transition. In addition, they can be produced as a result of v T decay at a relatively late stage. For 
the latter r d j Cay — r„ and for the former one should naively expect rj he ™ = (1/2) (4/3) = 2/3. If this is 
true the majoron mass should be too small to be interesting for the process of structure formation, since the 
characteristic scale of the structures to form first would be too large, as in the hot neutrino scenario. However, it 
is important to note that the majorons are rather strongly self-coupled particles, as they possess an interaction 
A J 4 whose coupling constant A could be as large as 10~ 2 . This interaction in the second order leads to the 
process of cannibalism, when collisions of four J's produce only two. The freeze-out of species whose number 
density was reduced by multiparticle collisions n — > 2 was considered in refs. ]l7| , |l8[ ] in a simplified approach. 

+A more detailed analysis of the nucleosynthesis bounds on u T masses and annihilation cross sections will be presented elsewhere 



Here we have numerically integrated the kinetic equation governing this process, under the usual assumptions 
that majorons are in kinetic equilibrium and described by Boltzman statistics. Then in complete analogy to 
two-body annihilation (see e.g. Il9[) we get 



where N is the combinatorial factor related to identical J's (we take N = 2), K\ is the McDonalds (Bessel) 
function, = nj/no, r e f = n e f/nQ, where n e } q is the equilibrium number density of majorons given by the 
Boltzman distribution with temperature Tj. 

It is a rather tedious job to calculate the cross section of the process 4J — > 2J but, since there is 
considerable freedom in the choice of the value of A which is not known anyhow, we make a simple estimate of 
the cross section tr(s) assuming that the amplitude of the process is a constant, A = X 2 /m 2 j. Since the photon 
temperature of is not equal to the temperature of the majorons, we need another equation which is provided 
by the covariant law of majoron energy density conservation. 

In our numerical estimates we have assumed the following parameter values: r v my[ = 1CP 2 , tjim = 10, 
T v — 5 yr., and flh 2 = 0.5. These numerical values are only illustrative, we use them as a reference set in 
order to show that the model can be self-consistent. With these numbers we get T\ » 630 eV, ~ 50 eV, and 
T2 « 0.45 eV. The frozen number density of majorons is found to be sufficiently small, rj w 0.2. We assume 
that the majoron mass is tok = 0.3. Correspondingly the characteristic size of the structures first formed can 
be two orders of magnitude smaller than in the neutrino universe, i.e. M = a few xl0 12 M©. 

Some bounds on the properties of MeV tau neutrinos and majorons can be found from astrophysics. 
In order to avoid excessive stellar cooling through majoron emission their mean free path inside a star (lj Vee ) 
should be smaller than the radius of the latter. For a supernova with temperature T k 10 MeV the latter is 
l*/" = \jan v = 10 10 (10- 5 /ff) 4 cm. For g > 10 we are safe. Another possible restriction arises from the 
observed neutrino flux from SN87A. There could be delayed neutrinos originated from the decay v T — > J -\- v 
p0| . But, for example, in the model of ref. Q the main v T decay channel is into so that this restriction is 
also harmless. 

By assumption the majoron arises from the spontaneous violation of the global U(l) lepton number 
symmetry which generates the neutrino masses. For v r masses in the MeV range the majoron has a sizeable 
diagonal Yukawa coupling to the v T of order g — 10 -4 . The majoron also has a much weaker nondiagonal 
coupling to muon neutrinos which induces the decay v r — > + J with lifetime of order of years. In order to 
avoid the supernova bound we assume the nondiagonal coupling to v e to be negligible. Moreover the diagonal 
couplings of the majoron to v e and are severely restricted in order to ensure that the majoron decay lifetime 
J — ► vv is larger than the age of the universe, i.e. g vu < 10~ 17 . Using the estimate m„ = g vll V and a typical 
value for the lepton number violation scale around 10 3 GeV we obtain that the v e and masses will be around 
10 -5 eV, i.e., just the range where the solar neutrino deficit will be explained by long wavelength or just-so 
v e to Un oscillations. 

Moreover, our unstable dark matter majoron will have tiny loop-induced decays to two photons J — > 77 
which will lead to the existence of an X-ray line at E 1 = mj/2 which could be detectable |^2[ . 

The self-coupling AJ 4 between majorons with A in the range from 10 -2 to 10 -4 leads to an interaction 
cross section a e i = \ 2 / (Q4irm 2 j ) which will maintain kinetic equilibrium between majorons till the present day. 





The corresponding majoron mean free path lf ree — (c e /«j) 1 with nj — n cosm = 10Vth 2 rn^ /cm 3 is about 100 
Kpc(mj<- / 'A|) where rax = mj/KeV and A3 = 10 3 A. Inside galaxies with nj — 3 x 10 4 cm _3 m^ 1 this mean free 
path is much smaller than the galactic size, I free = ^i m K / ^3) pc. 

The possibility that dark matter can have a relatively strong self-interaction, while being essentially 
decoupled from usual matter, has been considered earlier jl^, The behaviour of such matter is intermediate 
between that found in the CDM and HDM models. Structure formation with self-interacting dark matter 
particles has been recently criticized [ p2[ . Our scheme however is different from the simple original version and 
avoids at least some of these criticisms. The main objection of ref. is that it is impossible to generate 
enough small scale power to account for the damped Lya systems without producing too much power at cluster 
scales. Any model with decaying particles and early MD stage, such as ours, automatically takes care of this, 
since it has a suppressed power at cluster scales and simultaneously a large power on smaller scales due to rising 
fluctuations at the earlier MD stage. Another point raised in ref |^2| is that the galaxy merging should be 
quite different from what is expected in the usual collisionless dark matter case, and that the motion of galaxies 
with respect to the cosmic background would result in the stripping of the dark matter halos from galaxies. 
However, these points deserve further investigation. For example, the stripping of the halos may not take place 
if all the dark matter is clustered inside the halos and the cosmological background is absent. Present data 
show that the galactic halos extend up to two or more hundred kiloparsecs, and the 100% clustering at a large 
scale could account for Q w 1. Note also that the process of clustering of self- interacting matter differs from 
that of collisionless or baryonic matter because the former does not loose energy so easily by radiating photons. 

The self-interacting dark matter would form an isothermal matter distribution in the halos, and since 
the self-gravitating isothermal gas has a mass density which goes down with distance as 1/r 2 , it gives a natural 
explanation of the observed flat rotational curves. It was recently noted [^3) that the shape of the rotation 
curves inside or not far from the luminous center (for dwarves and also possibly for spirals) does not agree with 
the assumption of the collisionless dark matter. While this statement depends on how matter is distributed in 
galaxies, if confirmed, it may be a strong argument in favor of self-interacting dark matter. Another test of the 
model (suggested to us by S. White) is the shape of the halo: for self-interacting dark matter it should be more 
spherical than for the sterile one. 

Last but not least, the scenario described here can be implemented rather naturally in majoron extensions 
of the standard model of particle physics. As an example, supersymmetry with spontaneously broken R parity 
provides a consistent model for our scenario which incorporates all the necessary ingredients. 
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